function [ loglik_f ] = loglik_f_power( apost,apre, q_bar, qmin, betaq, sig, k)
%loglik_f Log likelihood of the asset density only
%   Detailed explanation goes here



q_up = exp(q_upperbar_proportional(k, q_bar ,[]));

F_pre   = log(f_density_power(apre, q_bar, qmin, betaq, sig, q_bar));
F_post  = log(f_density_power(apost, q_bar, qmin, betaq, sig, q_up));


loglik_f = [F_pre;F_post];

end


